Ultracold two-component Fermi gases with a magnetic field gradient near a Feshbach 
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We study theoretically the ultracold two-component Fermi gases when a gradient magnetic field 
is used to tune the scattering length between atoms. For ^Li at the narrow resonance Bq = 543.25 
G, it is shown that the gases would be in a coexistence of the regimes of BCS, BEC and unitarity 
limit with the present experimental technique. In the case of thermal and chemical equilibrium, 
we investigate the density distribution of the gases and show that a double peak of the density 
distribution can give us a clear evidence for the coexistence of BCS, BEC and unitarity limit. 
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With the remarkable development of the cooling tech- 
nique for two-component Fermi gases and Feshbach res- 
onance, the recent experiments j3i S 13 have finally re- 
alized molecular Bose-Einstein condensation (BEC) and 
the condensation of fermionic atom pairs on the side of 
attractive interaction has also been investigated by sev- 
eral experiments H H S H i, H IllJli- 

Together with 

intensive theoretical investigations [l2llldll4|.ll5l[T^ E^, 
Ifl, Jb 22, 23, 24, 25tEappll2ll2l IM M 

such as resonance superfluid jla, [l3 and universal be- 
havior for the gases with divergent scattering length 
[25L l2^ l2^ I28L 12a. I30I I , these experimental advances pro- 
vide us a quite unique system which will even contribute 
largely to our understanding on the mechanism of high- 
temperature superconductors. 

For two-component Fermi gases, when the magnetic 
field is tuned so that the energy of a quasibound molecu- 
lar state in a closed channel matches the total energy in 
an open channel, there is a magnetic-field Feshbach reso- 
nance [25I [3^ which can tune the scattering length from 
positive to negative over many orders of magnitude. On 
the side of repulsive interaction (BEC side), there exists 
molecule which is a short-range fermionic atom pairs. On 
the side of attractive interaction (BCS side), there would 
be a superfluid behavior due to the atomic Cooper pairs 
at sufficient low temperature. On resonance, the abso- 
lute value of the scattering length is much larger than 
the average distance between atoms and one expects a 
universal behavior for the system in the unitarity limit. 
In the present experiments on the BCS-BEC crossover, 
a uniform magnetic field is used to tune the scattering 
length through the magnetic-field Feshbach resonance. 
By tuning the uniform magnetic field near the resonant 
magnetic field, several experiments have investigated the 
BCS-BEC crossover such as condensate fraction 0, |3, 
collective excitation |3 , pairing gap ^3 1 and heat ca- 
pacity 01 . 

For ultracold two-component Fermi gases confined in 
an optical trap with axial symmetry along the z axis, in 
the present work, we will investigate the unique prop- 



erty of the system when the magnetic field to tune the 
scattering length has a gradient of a. For ^Li at the nar- 
row resonance Bq = 543.25 G, our researches show that 
the gases can be in a coexistence of the regimes of BCS, 
BEC and unitarity limit with appropriate parameters. 
For the gases in thermal and chemical equilibrium, when 
the pair size is much smaller than the cloud size of the 
system, the density distribution of the gases is calculated 
based on the local density approximation. 

For Fermi gases with an equal incoherent mixture of 
the internal states |1) and |2), near a Feshbach reso- 
nant magnetic field Bq, the scattering length is a [B) — 
AaQ (1 — w/ {B — Bq)) with gq being the Bohr radius. 
For the magnetic field in a;— direction with gradient a 
along z— direction, we have B — {Bq + az) Cx- fn this 
case, the scattering length becomes z— dependent and 
takes the following form: 



a (z) — Auq (1 
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ff the magnetic field gradient a is positive, when z < 
and a (z) n^/'^ < 1 (n is the total density distribution 
of the ultracold gases), there is repulsive interaction be- 
tween fermionic atoms which corresponds the BEC side, 
while when z > and |a(z)|n^/'^ < 1 there is attrac- 
tive interaction between fermionic atoms which corre- 
sponds to the BCS side, fn the regime determined by 
|a(z)|n^/'^ > 1, the ultracold gases are in the unitarity 
limit where the scattering length can be regarded to be 
divergent and will not appear in the final result of a phys- 
ical quantity such as chemical potential. Thus, by using 
the gradient magnetic field and appropriate parameters, 
there would be a coexistence of the gases in the form of 
BCS, BEC and unitarity limit. 

For the optical trap with angular frequencies being 
ujr{— = ^y) and ujz in the radial and axial direc- 
tions, when the magnetic field with gradient a is applied 
to the ultracold gases, the overall external potential for 



2 



the fermionic atom in the internal state \i) is 



2 / 2 I 2\ , 2 1 



(2) 

where i = 1,2. Due to the presence of the magnetic 
field, for different internal state |z), there is a differ- 
ent shift to the optical trap which is determined by 
Az|i^ — a^i^^^^ jmuJ^ with being the magnetic mo- 

ment of the internal state In the above expression, 

AV|i) — —a^ ('^iT)"^) I'^rnujl- One can also get the over- 
all external potential for a molecule 



(3) 



where /i™"^ is the magnetic moment of the molecule. 

The regime of the gases in the unitarity limit can be 
roughly determined by \a {z) \ v}!'^ > 1. Thus, the length 
scale for the regime in the unitarity limit can be esti- 
mated as 211'^'^ with l^'^'^ — Awn^^'^ao/a. Because there 
is a different shift for the external potential of the inter- 
nal state |1) and |2), the difi'erence of the shift should be 
much smaller than 211^'^ so that one can omit safely the 
effect of the shift Az\iy After a simple calculation, this 
condition |Az|i^ — Az|2)| << 211*^^ means that 
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the average distance between atoms |3J| . For the atomic 
Cooper pairs on the BCS side, one can estimate the pair 
size based on the BCS theory. We have verified that for 
'the parameters used in the present work, the pair size is 
much smaller than the cloud size. Thus, one can safely 
use the well-known local density approximation to calcu- 
late the density distribution of the system. 

To calculate the density distribution of the system, we 
firstly investigate the chemical potential of the system. In 
the regime of molecular BEC (a (z) -n}/^ < 1) with den- 
sity distribution n^^'^ and molecular scattering length 
o-m = 0.6a !3^. the chemical potential of the molecular 
BEC is given by 
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where Afi — — ^ . The last term AfiBo in the 

above equation is due to the energy of the bound state in 
the closed channel when there is no magnetic field. After 
a simple calculation, one can get the following expression 
for the chemical potential 



r)_fc2 „BEC 
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where e,„ — — with ~ f? jmcP' being the binding 
energy of the molecule. V^j^; is given by 



For the typical parameters that uiz/2tt ^ 100 Hz, n ^ 
IQi^ cm~^ and V^'^'^ ~ 10 /Ltm, the condition given by Eq. 
Q requests that Aw is smaller than 10 G. This means 
that very narrow resonance for an element is an appro- 
priate choice to investigate the coexistence of the regimes 
of BCS, BEC and unitarity limit. A careful investigation 
shows that the mixture of |1) = \F = 1/2, mp = 1/2) 
and |2) = |F = 1/2, mp — —1/2) for ^Li at the narrow 
resonance located at Bq — 543.25 G ^3 '^^^ satisfy our 
request. For this narrow resonance. Aw is estimated as 6 
G [33j. Thus, in the present work, we will consider two- 
component Fermi gases of ^Li at the narrow resonance 
Bo = 543.25 G. 

For two-component Fermi gases of ^Li confined in the 
optical trap, we will investigate the density distribution 
of the gases by using the condition that the system is in 
thermal and chemical equilibrium. Our researches show 
that with appropriate choice of parameters and element 
for the system, the coexistence of the gases in the form 
of BCS, BEC and unitarity limit can be clearly shown 
through the density distribution. In the present exper- 
iments, the gases can be cooled far below the critical 
temperature of molecular BEC and Fermi temperature 
of the Fermi gases. Thus, we will consider the system 
at zero temperature to give a clear presentation. For 
the molecules on the BEC side, the pair size can be re- 
garded as the scattering length. For the atomic pairs 
in the unitarity limit, the pair size can be estimated as 



-E,M";'Bo + Ayw, (7) 

where Az^o; = a (^E^^""^^ /2mujl and AV,noi = 

-a^ (E,Mirr)'/4mc.2. 

In this regime, the chemical potential of the Fermi gas 
in the internal state \i) takes the form 
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where n^^f)^ is the density distribution of the Fermi gas 



i). In the above expression gaa = 
: m/2 being the reduced mass. In 



in the internal state 
2Trh? a /rrir with rrir 
addition, gam — 0-Qgaa [2^1 which is obtained based on 
the atom- molecule scattering length Oam = l-2a |35l |. 

In the regime of the unitarity limit where the absolute 
value of the scattering length is much larger than the av- 
erage distance between atoms, as pointed out in [33|, we 
assume that the gases is in the mixture of Fermi gases and 
dimeric gas in chemical equilibrium. Based on the local 
density approximation at zero temperature, the chemical 
potential of the Fermi gas with density distribution n^^ 



takes the form 
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where j3i is firstly measured in [25j and /9i = —0.56 based 
on a quantum Monte Carlo calculation |27j • Omitting the 
binding energy of the dimer in the unitarity limit, based 
on the dimensionality analysis, the chemical potential of 
the dimeric gas with density distribution n^^ can be as- 
sumed as 



2 X 2m 



2/3 



Vd. 



(10) 



For the dimeric gas, Vd takes the same form as V^^i 
for the molecular gas. Based on the condition of the 



chemical equilibrium Si/i^i^^ = Md^i one can get the ra- 



tio = [4(l + /3i)/(l + 



From the experimental result n^i^ /vH^ 



3/2 



on resonance. 
4 in 0, (32 is 

estimated as —0.3. 

In the regime of BCS, the chemical potential can be 
approximated as 
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where n^R^ 
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is the density distribution of the fermionic 
atoms on the BCS side. 

In the case of thermal and chemical equilibrium for 
the system, the minimum of the Gibbs free energy means 
that there is following important relation for the chemical 
potential: 
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(12) 

To illustrate clearly the density distribution and evidence 
for an experiment to show the coexistence of the regime 
of BCS, BEC and unitarity limit, we use the parameters 
tj^/27r = 120 Hz and vJjyl^ = 0.2 x 10^^ cm"^ at the center 
z — Q. In addition, the gradient of the magnetic field is 
chosen as a 
|Az,,)|/2/-^ 



21 G/m. Based on these parameters. 



« 1, |AV^,)|/^ « 1, Az 

mol 



and 



AVmoi/fJ- = 0. Thus, one can safely omit A2|j^ 
Azmoi and AVmoi- 

From the equilibrium condition for the chemical poten- 
tial given by Eq. (|12|l . one gets the following expression 
for the density distribution: 
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(13) 
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FIG. 1: Shown is the density distribution of the gases when a 
gradient magnetic field of 21 G/m is used to tune the atomic 
scattering length. After the molecules in the regime of BEC, 
the dimers in the regime of unitarity limit, and the atomic 
Cooper pairs in the regime of BCS are converted into deeply 
bound molecules by decreasing non-adiabatically the mag- 
netic field, the thick dashed line illustrates the density dis- 
tribution of the deeply bound molecules. We see that there 
is a double peak for the density distribution of the deeply 
bound molecules which can give us a clear evidence for an 
experiment to make the system in the coexistence of BEC, 
BCS and unitarity limit. The ratio of the peak density to 
valley density is estimated as 1.4. In this figure, the density 



distribution is in unit of n 



fW 



(^ 



0), while the coordinate z 

is in unit of Rz = \/ ^^f f /mujl. In the inset of this figure, 
shown is the ratio R of the peak density due to molecular 
BEC and ultracold gases in the unitarity limit for different 
magnetic field gradient. 
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where Vext takes the form 



Vext = \rnLol {x^ + y^) + ^mujlz^. (18) 



To give a concise presentation, we have introduced 
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Using the parameters in this paper, shown in Fig.l 
is the density distribution of the gases where the den- 
sity is in unit of n^f^i^ (-^^ — 0), while the coordinate z is 
in unit of Rz — -^//ie// /mii;2 . The arrow shows the lo- 
cation ±0.7i?2 which is determined by |a (z)] n"'^/'^ = 1. 
The solid line in Fig.l shows the density distribution 
of molecular BEC nf^^^ (z < -0.7), while the dot- 
dashed line in this figure shows the density distribution 
n^^ (—0.7 < z < 0.7) of the dimers in the unitarity 
limit. The dotted line shows the density distribution 
of nfy^f {z > 0.7). Analogously to the experiments in 
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, we consider the process that the magnetic field is 
decreased non-adiabaticaUy so that the fermionic atom 
pairs (i.e. molecules on the BEC side, dimers in the uni- 
tarity limit, and atomic Cooper pairs on the BCS side 
before the decreasing of the magnetic field) are converted 
into deeply bound molecules. After this non-adiabatical 
process, from the density distribution of molecular BEC 
(solid line) and dimers in the unitarity limit (dot-dashed 
line) , we see that the coexistence of the gases in the form 
of BCS, BEC and unitarity limit can be clearly shown 
through the double peaks in the density distribution of 
the deeply bound molecules. The thick dashed line shows 
a schematic interpolation of the density distribution of 
the deeply bound molecules. The evidence for the coex- 
istence of the regimes of BCS, BEC, and unitarity limit 
can be also shown through the non-symmetric density 
distribution of the deeply bound molecules about z = 
after the non-adiabatic process. The inset in the figure 
shows the ratio R between the peak density of the regime 
of BEC and unitarity limit for different magnetic field 



gradient. Through the density distribution for different 
magnetic field gradient, our researches show that for 18 
G/m< a < 22 G/m, there is obvious double peak density 
distribution. 

In summary, we show that for ^Li at the narrow res- 
onance Bq — 543.25 G, by using a gradient magnetic 
field to change the scattering length, one can make the 
gases become the coexistence of the regimes of BCS, BEC 
and unitarity limit. In the case of thermal and chemical 
equilibrium, it is shown that with appropriate parame- 
ters there is a double peak in the density distribution of 
the deeply bound molecules after the non-adiabatic de- 
creasing of the magnetic field. This can give us a clear 
evidence for the coexistence of the regimes of BCS, BEC 
and unitarity limit in an experiment. 
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